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, Abstract. We analyze free-fermion conditions on vertex models. We show -by 
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mann weights, and propose a general scheme for such a process in two and more 
dimensions. 
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1 Introduction 



We consider vertex models in D dimensions. As is standard, the bonds of 
the lattice carry variables taking q values (colors). The model is determined 
by attributing Boltzmann weights to the various possible bond configurations 
around a vertex These homogeneous weights are arranged in a matrix, 
which we denote by R. The size and form of the matrix R vary according to the 
number of colors, and the coordination number of the lattice. Typical examples 
we will consider are the 2D square lattice, the 2D triangular lattice, and the 3D 
cubic lattice: 




2D square 2D triangular 3D cubic 



If the number of colors q is 2, and we will restrict ourselves to this case, then 
the i?-matrices are of sizes 4 x 4, 8 x 8, and 8x8 respectively. The difference 
between 2D triangular and 3D cubic for example does not show in the size nor 
the form of the matrix. It will, however, appear in the operations we define on 
the matrices. 

We shall use a number of elementary transformations acting on the matrices. 
These transformations come from the inversion relations and the geometrical 
symmetries of the lattice, in the framework of integrability ^ ||, [| M, and 
beyond integrability H . They generically form an infinite group r/ Qttice [Q, |[ . 

The groups Ti at uce have a finite number of involutive generators. The first 
one, denoted J, is non-linear and does not depend on the lattice: it is the matrix 
inversion up to a factor. The other generators act linearly on R, actually by 
permutations of the entries, and represent the geometrical symmetries of the 
lattice. 

For the square lattice, we have two linear transformations, the partial trans- 
positions ti and t r ||: 

{tiR)%=Rt, {t r Rfl=R% i,j,u,v = l..q (1) 

The product Ut r is the matrix transposition (7 stand for 'left' and r stands for 
'right' in the standard tensor product structure of R). 

For the triangular lattice, we have three linear transformations 77, T m , T r : 

(T,RYi k — R" fe j I U\ijk _ nwvu r RY3 k — R-?™ (9\ 
\'l ll )uvw ~ ll iwv \' mJ-lJuvw ~ ll kji > I'r^iiimii — lx vuk' \* I 
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Finally, for the cubic lattice, we have three linear transformations ti, t m , 
tr §: 

and the product ti t m t r of the three partial transpositions is the matrix trans- 
position. 

All the generators are involutions. All products tjj, t m I and t r I are of 
infinite order when acting on a generic matrix, as are T(I and r r I. On the 
contrary r m and / commute and (r m /) 2 = 1. 

It is straightforward to check that: 

til ti = ti I ti and t r I t r = r r I r r (4) 

so that essentially Y tr iang appears as a subgroup of Y cu uc, up to finite factors. 

It is important to keep in mind what the "size" of the groups Y are. All 
three Y square , Y trtang , and T cubic are infinite, but T square has one infinite order 
generator, Ttriang has two, and T cu bi C has three. The last two groups are thus 
hyperbolic groups , and studying the triangular lattice can be a good test-case 
for the more involved tridimensional cubic lattice. 

The groups Y are the building pieces of the group of automorphisms of 
the Yang-Baxter equations and their higher dimensional generalizations, and 
solve the so-called "baxterization problem" || . These equations form overde- 
termined systems of multilinear equations, of which the possible solutions are 
parametrized by algebraic varieties S. The overdetermination increases very 
rapidly with the dimension of the lattice. At the same time, the size of Y also 
explodes. When looking at solutions of the Yang-Baxter equations and their 
generalizations to higher dimensional lattices, one faces a conflict between hav- 
ing a more and more overdetermined system and a larger and larger group of 
automorphisms for the set of solutions. We will show how this conflict is re- 
solved in some 2D and known 3D solutions by a degeneration of the effective 
realization of the group T, which becomes finite. 

The content of this letter is the description of a mechanism for such a de- 
generation, obtained by the linearization of specific elements of Y. 

We first show how the free-fermion condition on the asymmetric eight vertex 
model Q falls into this scheme. 

We then describe the group T tr iang for the 32- vertex model on the triangular 
lattice. We show that the free-fermion conditions given in (To[ |ll[ amount to 
linearizing the inverse I and make the realization of Y tr iang finite. 

We finally write and discuss similar conditions for the 32-vertex model on 
the cubic (3D) lattice, by analyzing solutions of the tetrahedron equations jl2]. 

[IB- 

One of the results we obtain is that free-fermion conditions should always 
appear as quadratic conditions, whatever the size and form of the matrix R is, 
and in particular whatever the dimension and geometry of the lattice are. 
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There already exists an important literature about free-fermion models. We 
may refer to |l5|, [l6|, [l7|, where an exploration of the use of grassmannian 
variables, both for the construction and the resolution of the models, can be 
found. This work also motivated the interesting 3D construction of p8| . 

Our approach is based on a direct study of the matrix of Boltzmann weights, 
concentrating on the action of the symmetry group T, and provides another view 
on this class of models. 



2 Some notations 

At this point it is useful to introduce some notations we will use in the sequel. 

We will denote the equality of two matrices R and R' up to an overall factor 
by R ~ R'. 

We always denote by t the full matrix transposition. 

We will use various gauge transformations (weak graph dualities) |1| , that 
is to say the conjugation by invertible matrices which are tensor products, also 
defined up to overall factors, i.e. transformations of the type 

R — > gi 1 <8> ® ■ ■ ■ ® 9k 1 ■ R ■ 9i ® 92 ® ■ ■ • ® gu ■ 
Define the matrices 

CTo= (o ?)' CTl = (i o)' CT2= (i o 1 )' 173 = (o -i) (5) 

and the matrices <r ai a 2 ...a k of size 2 k x 2 k (k will be 2 for the square lattice, 3 
for the triangular and cubic lattice, and so on) by: 

& ai a 2 ...a k = (J ai ® <?a 2 ® ■ ■ ■ a a k ■ 

We denote by S ai a 2 ...a fc the conjugation by a ai a 2 ...a k ■ 

Clearly both t and I commute with all S aia2 ... afc , up to an irrelevant sign. 
Moreover, from the fact that o a &b = ± &b &a, Va,6 = 0, 1,2,3, the gauge 
transformations £ satisfy 

S ai o 2 ...a fc ^b 1 b 2 ...b k = i ^bib 2 :-b k ^- , axa 2 — a k i 

meaning that they commute up to a factor. 

Particular gauge transformations of interest are 

= S33...3, 

and some transformations acting just by changes of sign of some of the entries, 
and denoted e a (a = I, m, r, ...). 

If fc = 2 (square lattice), then a = I or r, and 

e i = S30, e r = £03. 
If k = 3 (triangular and cubic lattice) : 

e l = £300: e m = £o30i £ r = £(303 ■ 
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3 Free-fermion asymmetric eight-vertex model 

The matrix R of the asymmetric eight- vertex model [^0| is of the form 

(a d'\ 

R= c b> (6) 
\d a'/ 

Notice that this form is the most general matrix satisfying it R = R. 
The free-fermion condition || (see also j2lj) is 

aa' - dd' + W - cd = (7) 

A matrix of the form (^J) may be brought, by similarity transformations, to a 
block-diagonal form 

*-(* »), with *-(; *) and «.-(* ^ 

If one denotes by 8\ — aa' — dd' and by 62 = bV — cd the determinants of 
the two blocks then the matrix inverse / written polynomially (namely R — > 
det(R) ■ R^ 1 ) just reads 

a — * a' ■ 82, a! — ► a ■ 82, d — > — d -82, d 1 — > — d' • ^2) 
& — ► 6' ■ Si, —* b • 81, c — > — c • 5i, c' — » — c' • Ji. 

The condition (0) may be written as pg(R) = with the notations of ||, and is 
consequently left invariant by r sguare . It is straightforward to see that condition 
(0) is 81 = — 82 and has the effect of linearizing I into 

a — ► a', a' — ► a, d — » — d, d' — > — d', 
6 -> -6', 6' -> -6, c -> c, d -> c'. 

The group T is then realized by permutations of the entries, mixed with changes 
of signs, and its orbits are thus finite. The commutators of partial transposi- 
tions and inversion, in the sense of group theory, i.e: tiltY 1 !^ 1 = (U I) 2 and 
trlt^ 1 !^ 1 = (t r I) 2 reduce to a change of sign of the non-diagonal entries of R. 
These commutators are typical infinite order elements of T, when acting on a 
generic matrix, and their degeneration is a key to the finiteness of the realization 

of r. 

If we introduce the grading gr 



gr(R) 





-b -d 

-c -b' 
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which operates by changing the sign of the entries of only one of the two blocks, 
say i?2, then, for any R satisfying (R), the action of the inverse reduces to 



I(R) ~ t Ei 2 gr{R) 



(8) 



where t is matrix transposition. In other words we have defined, on all matrices 
satisfying n(R) = R, a linear operator 



l sq = t E 12 gr 
such that the free-fermion condition (0) reads 

I(R) ~ 

or equivalently 

i? • l sq (R) — unit matrix 
The linear transformation Z S(? satisfies a number of relations: 



id, 



J, r 



(9) 
(10) 

(11) 



Such relations ensure that the orbit of R under T is finite, as is readily checked, 
and specify the changes of signs to which (ti I) 2 and (t r I) 2 reduce. Notice that 
the definition of l sq is not unique. 



4 Free-fermion conditions for the 32- vertex model 
on the triangular lattice 

We consider the free-fermion conditions for the 32-vertex model on a triangular 
lattice, and use the notations of Q: 



R 



fo 








$23 







fl2 








f36 


fm 





fl6 













f35 


fz5 





flB 








he 


f56 













fn 


/s^ 








fs4 


fs4 















he 








fl5 





jW 


fs5 





fi5 








fl6 





/S6 


JS6 








fis 


fl3 





J23 








fo 



(12) 



This matrix may be brought, by a permutations of lines and columns, into a 
block diagonal form: 

i?i 
R 2 



i? = 



with 



(13) 
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Ri = 



Ri = 



fi4 fa4 {24 fee 

fl6 fad fae fj5 

fl5 }35 fa-5 f{6 

fas }l2 flS fa 



(14) 



fa fl3 fa 2 fj>3 

fae fas fss flS 

fa5 fad fs6 fl6 
fad fa4 fs4 fl4 

The inverse / written polynomially is now a transformation of degree 7. If one 
introduces the two determinants Ai = det(Ri) and A2 = cfei(i?2), then each 
term in the expression of I(R) is a product of a degree three minor, taken within 
a block, times the determinant of the other block. 

Denoting /12 = fa^e and so on, the free-fermion conditions of 1 1 , [ll] arc: 



fafijkl fijfkl fikfjl filfjhj 

fofo = /12/12 — /13/13 + /14/14 



\/i,j,k,l = 1,. 

/l5/l5 + /l6/l6 



,6 



(15) 
(16) 



What is remarkable is that, not only the rational variety V defined by fllqjlq ), 
is globally invariant by Ttriang, but again the realization of T on this variety is 
finite. This comes from the degeneration of / into a mixture of changes of signs 
and permutations of the entries, as was the case in the previous section. 

When relations (|l5|, [ll]) are satisfied, the action of / simplifies to 



with 



I(R) ~ l tr (R) 



hr — t S121 ~ t Em e ri 



(17) 



(18) 



Since we have the prejudice that all free-fermion conditions should be invari- 



ant under T tr iang, as (15, |16|) are, one should complement (17) with 

I n(R) - hr n(R) 

lT r {R) ~ hrT r (R) 

We may list some useful relations: 

77 2-labc — ^acb ^~rn ^abc — ^c6a 7"mi T r Yj aoc — Yj^ ac T r 

t r a = r a t V a, b, c = 0, 1, 2, 3, V a = I, m, r 

(ja Tf3 r 7 ) 2 = id Va, /?, 7 — l,m,r 

TlT m T r I = I TlT m T r , TlT r T t I = I TlT r T t 

The linear transformation L r satisfies in addition: 



72 

Hr 



hr tl hrU — £f, 



— ZG?, hr t — ^ hr 

hr tm hr^ra — %d^ h r t r hr^r 



(19) 
(20) 



(21) 
(22) 
(23) 
(24) 



(25) 
(26) 



Using relations ( ^l| ) to (p6|), it is possible to show that the completed system 
( 171 pi p0| ) is left invariant by the action of the group T tr iang- 
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The system is also invariant under the gauge transformations leaving the 
form ([12]) stable. (Hint: The gauge transformations leaving the form (|l^) stable 
satisfy 7 E a bc 7 — Sabc when a, b = 1, 2). 

Moreover the generic 32- vertex invertible solutions of the completed system 
(0 [l9|, |0|) satisfy the free-fermion conditions (p|, 

It is clear from (jl7|, [l9[ |2(]) an d (f2l|) t o ( |26| ) that the realization of the group 
T triang is finite, when conditions (|l5|, |16[) are fulfilled. 

One should also notice that any linearization condition of the type of ( [l7| ) 
is a set of quadratic conditions, whatever the size of the matrix is. Indeed they 
mean that the matrix product of R with some linear transformed C(R) of R is 
proportional to the unit matrix, i.e: 

R ■ C(R) ~ unit matrix (27) 

and this is a set of quadratic conditions. 

Remark: The invertible solutions of ( |l7j ) form a group for the ordinary 
matrix product, since I ■ lt r is an automorphism of the group of invertible ma- 
trices of the form (jT|), i.e. I{l tr (Ri • #2)) = I{kr{R\)) • I{kr{Ri))- The extra 
conditions added when completing the system break this in such a way that the 
ordinary matrix product of three solutions is another solution. In other words, 
if Ri, i?2, R3 £ V, then Ri ■ R 2 ■ R 3 £ V, while i?i • R 2 <£ V. This was actually 
already the case for solutions of (Q), but the mechanism is more subtle here as 
conditions (JTJ, [if]) imply Ai = +A 2 . 



5 32-vertex model on the cubic lattice 



We now turn to a solution of the tetrahedron equations [jl4| , 
of the form 



R 



( 


d 








—a 





-b 


c 












w 


a; 





y 

















X 


w 





z 








y 






—a 








d 





c 


-6 












-2/ 


z 





ic 








— X 






b 








c 





d 


a 









c 













a 


d 







V 





z 


-y 





—x 








w 


) 



Let i? be 



(28) 



The form of ( ]2q ) is stable under the group r cu ;,ic, and it is natural to look for 
invariants of T in the space of parameters {a, 6, c, c?, a;, y, z, w} 24 . There exist 
five algebraically independent quadratic polynomials in the entries, transforming 

1 Notice that the form (|2?|) is not stable by the circular permutation of the three spaces 
{l,m,r}. 
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covariantly, and with the same covariance factors under all generators of T cu bi C - 
They are: 

ax, by, cz, dw, and Q = a 2 + c 2 — d 2 — y 2 — b 2 + x 2 — w 2 + z 2 . 
We thus have four algebraically independent invariants of T cu bic, say for example 

ax by cz Q 

Xi = -1- j X2 = -r- , X3 = -j— , and xo = -5- ■ 

dw aw dw dw 

A complete analysis shows that there is no other algebraically independent in- 
variant of r 'cubic- A numerical and graphical study p5[j , shows how "big" the 
realization of T is for generic values of the above invariants. 



These invariants are completely specified in the solution 14 , for which 



Xi = X2 = X3 = 1 (29) 
Xo = (30) 

Out of the four invariants, xo plays a special role. If xo = 0, then the action of 
I linearizes quite in the same way as in the previous cases. Notice that, strictly 
speaking, condition ( p0| ) is not so much an assignment of value to the invariant 
Xo but rather a vanishing condition for the covariant quantity Q. Recall that 
assigning a definite value to an invariant object is meaningful whatever this 
value is. On the contrary covariant objects cannot be assigned a value unless 
this value is zero. 

When Q = 0, one gets 

I{R) ~ l c (R) (31) 
The linear transformation l c may be written 

l c = t £ 03 o 9 r i 

where t is transposition and gr is a grading changing the sign of the entries of 
R belonging to the same block, say {x, y, z, w}. Notice that the definition of l c 
is not unique, due to the very specific form of (p8|). Notice also that Q = is 
one of two quadratic conditions ensuring the equality of the determinants of the 
two blocks of R (see (|l3|)). The other one is not stable under T. 
The linear transformation l c satisfies 

I 2 =id, l c t = t l c , (32) 
lc t a l c ta = f a , \/a = I, m, r (33) 

Any matrix of the form (|2^) with Q = obeys 

I(R) ~ l c (R) (34) 
It a (R) ~ l c t a (R) Va = l,m,r (35) 
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Using (32 3^), it is straightforward to show that the complete system (|34|,[35|) is 
invariant under T cu uc and that the orbit of R is finite.^ 

The study of the additional conditions (^9|) would take us beyond the scope 
of this letter, but we may make a few remarks. 

The first remark is that since among conditions (|2^,|3(]) , only (|3^) has to do 
with the fmiteness of the realization of T, ( p9|) may have nothing to do with 
free-fermion conditions. They are additional constraints making the resolution 
of the tetrahedron equations possible, and this may be understood as follows. 

The tetrahedron equations are in essence a compatibility condition for the 
existence of non-trivial solutions of the "propagation properties" p^] (alias 
"Zamolodchikov algebra" |2j| , alias "vacuum curves" |2^, Q , alias "pre-Bethe 
Ansatz" equations |3[ j30f): 

What conditions ( 29 3C|) e nsure is the existence, for fixed R, of a one-parameter 



family of solutions of (|3q). In the case we consider here, the family happens to 
be parametrized by a curve of genus larger than one. 

By eliminating {q, q 1 , r,r'} (resp. {p,p',r, r'} or {p,p',q, q'}) from ([36]), one 
gets conditions relating p,p', (resp. q, q' and r, r'). Such relations are generi- 
cally of degree 8 (biquartics). One effect of ( |2^ , |30| ) is that they all reduce to 
asymmetric biquadratic relations, defining three genus one curves of the form 

xV - i + (y 2 - x 2 ) k xv = (37) 



with 



be (d 2 - a 2 ) _ ac jb 2 + d 2 ) _ ab {c 2 - d 2 ) 



ad(b 2 -7~Y qq '~ bd (a' + c 2 )' ' urr ' ~ cd(a 2 -a 2 )' 



These three elliptic curves have different (algebraically independent) moduli. 
Their asymmetric character may be an obstacle to the use of ( |36| ) in the con- 
struction of the Bethe Ansatz states [J since the composition of relations of type 
([57]) reproduces the same type of relations, but alters the value of n by: 



2 \ k 



(38) 



Exceptional values of k (±1,oo), yielding a rationalization of (^), are fixed 
points of d38|). For these exceptional values, in particular k ~ 00, obtained with 
d = 0, the construction of a 3D Bethe Ansatz may be envisaged. 

?This is also the case for the bidiagonal solution of the "constant" tetrahedron equations 
^R.J. Baxter, private communication. 



6 Conclusion 



We have shown, through specific examples, how free-fermion conditions turn 
into degeneration conditions of our groups T: the generically non-linear (ratio- 
nal) infinite realization of T becomes a linear finite group. 

We believe this is a characteristic feature of 2D free-fermion models. 

We have shown that the known vertex solution of the tetrahedron equa- 
tions does have such a feature. An appealing issue is to decide whether or not 
such a statement can be made about other 3D and higher dimensional mod- 
els. Of course the full answer will come from linking directly the phenomenon 
we describe with explicit calculus using grassmannian variables. The particu- 
larly simple form of the conditions (combinations of products of entries with 
plus and minus signs), and the linearization process of the inverse should stem 
from elementary properties of exponentials of quadratic forms in anticommuting 
variables. 

Producing new solutions of the tetrahedron equations is another challenging 
problem. What could be done is to look for forms of the matrices R enjoying 
the linearization property we have described. This is a rather simple way to 
produce "reasonable" Ansatze for R. 

The next step would then be to study the so-called propagation properties 
(see above) rather than confronting directly the tetrahedron equations them- 
selves. 

Indeed these simpler equations, because they govern the construction of 
Bethe Ansatz states -a basic in the field-, underpin 2D, 3D, and higher dimen- 
sional integrability. 
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